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Phononic computing - the use of (typically thermal) vibrations for information processing - is a nascent 
technology; its capabilities are still being discovered. We analyze an alternative form of phononic computing 
inspired by optical, rather than electronic, computing. Using the acoustic Faraday effect, we design a phonon 
gyrator and thereby a means of performing computation through the manipulation of polarization in transverse 
phonon currents. Moreover, we establish that our gyrators act as generalized transistors and can construct 
digital logic gates. Exploiting the wave nature of phonons and the similarity of our logic gates, we demonstrate 
parallel computation within a single circuit, an effect presently unique to phonons. Finally, a generic method of 
designing these parallel circuits is introduced and used to analyze the feasibility of magneto-acoustic materials 
in realizing these circuits. 



PACS numbers: 85.70.Ec, 63.20.kk, 72.55.+S, 43.25.-x 

Computers are unquestionably useful, thanks in no small 
part to their silicon-based electronic architecture. After all, 
the size and expense of vacuum tube computers rendered per- 
sonal computing impossible. Today, there are many proposed 
alternatives to silicon computers. Quantum computing, for ex- 
ample, promises to perform calculations more efficiently than 
any classical computer could [1]. Moreover, biological and 
fuzzy logic computers outperform classical computers in spe- 
cific applications but are unlikely to supplant traditional com- 
puters [2, 3]. The creation of phonon (vibrational) equivalents 
of electronic computer elements has recently been explored. 
Phononic (often called thermal) computing's typical applica- 
tion is to use heat from electronic computers, where waste 
heat is cheap, plentiful, and wasted [4, 5]. 

The development of further applications for phonon com- 
puters has been stymied by the difficulty in controlling phonon 
propagation. This is aggravated by reliance upon analo- 
gies between temperature and voltage, thus taking inspiration 
from electronics [4-8]. Such focus has yielded reliance upon 
nanostructured materials [9, 10] and ID, nonlinear chains [4- 
6, 1 1] to construct circuit elements. While this has resulted in 
difficulties with constructing phonon computers, electronic- 
inspired computation has also impacted how information is 
encoded. Prior efforts encoded information in temperature 
differences. But when phonons alone are information carriers, 
one could take inspiration from optical computing, encoding 
information in the polarization of the phonon current. Log- 
ical 1 (0) would be vertically (horizontally) polarized trans- 
verse phonons. Logic operations in this computer therefore 
correspond to modifying the polarization of phonon currents. 
Kittel first realized the possibility of controlling phonon po- 
larization. He discovered the acoustic Faraday effect (AFE, 
or acoustic Faraday rotation, it rotates the angle of polariza- 
tion) for ultrasonic waves in ferromagnets subject to a uniform 
magnetic field [12]. Interest in this application seems to have 
ended there. Instead, the focus on magneto-acoustic (MA, 
also called magneto-elastic) effects (the influence of a mag- 
netic field on phonons) has been to discover new mechanisms 
for MA effects and how they might aid in characterizing ma- 
terials [13-15] or controlling magnons [16]. 

Here we consider the construction of logic from MA ma- 
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FIG. 1: (COLOR ONLINE) Comparing electronic and phononic 
computing. Electronic computing uses a stream of electrons through 
a solid-state transistor to control the magnitude of the voltage at the 
drain. Conversely, our implementation of phonon computing has 
transverse phonons undergo an acoustic Faraday rotation, controlling 
the polarization of the transmitted signal. Note the use two incident 
phonon currents with different polarizations and frequencies, each 
rotated by a different angle while traversing the gyrator. 



terials beginning with the fundamentals - sources, wires, and 
operators. Attention will be paid to gyrators, where we find 
that, in addition to isolators (the optical analog of diodes), gy- 
rators can create transistors. Logic gates are constructed from 
these transistors. We analyze how individual gates could per- 
form a new form of parallel computing (which we term rain- 
bow parallelization, RP) (Figure 1). 

It is helpful to consider where a signal comes from and 
how it reaches computing elements. Waste heat would be 
the ideal source of cheap, convenient phonons, but it lacks 



coherence, intensity, and polarization. At the other extreme, 
the phonon laser [17] has good beam characteristics but is in 
the early stages of development. For MA, acoustic transduc- 
ers are a promising compromise. Piezoelectrics (e.g. AC- 
cut quartz) convert electronic signals into polarized phonons 
currents, with polarization depending on the crystallographic 
plane at the interface [18]. At higher frequencies (GHz to 
THz, instead of MHz to GHz), femtosecond lasers can be 
transduced into picosecond phonon pulses via thin metallic 
films (e.g. polycrystalline Al on Zn) [19, 20]. Once phonons 
have been produced, they can be steered via wires. Any ma- 
terial with good lattice thermal conductivity would do, as this 
implies that phonons propagate quickly and retain polariza- 
tion over long distances (long phonon mean free path). 

While the previous elements are similar to their electrical 
counterparts, operators are not. In digital electronics, oper- 
ators (e.g. transistors) control voltage levels (logic states) 
by blocking or admitting electrical current; whereas opera- 
tors in polarization-based computing always admit currents 
(states are modified by switching polarizations). Given our 
basis, a change of state is accomplished by having a gyrator 
rotate the polarization by ir/2. Gyrators can be constructed 
from any material possessing the AFE (e.g. MA) (Figure 
1). In the AFE, a longitudinal magnetic field (B\\k, k is the 
phonon wave-vector) induces circular birefringence (two cir- 
cularly polarized modes travel at different speeds). Every lin- 
early polarized mode is composed of a superposition of cir- 
cularly polarized modes, which reach the end of the gyrator 
at different times. This phase difference between the circular 
modes yields a linearly polarized current with a new polariza- 
tion. The magnetic field strength determines the magnitude of 
the gyration [21]. There are many MA materials known [13], 
although yttrium iron garnet (YIG) is a good candidate here 
as the AFE has been observed at room temperature [21]. 

Currently, when gyrators are used, they are paired with a 
fixed magnet (often 0.01 -10T) to give a predefined gyration. 
This design has some uses in computing (e.g. when creat- 
ing isolators), but gyrators have other uses. We propose that 
a phonon current could control the strength of the magnetic 
field, allowing the gyrators to operate in a completely new 
way. One control mechanism is to transduce a second current, 
rectifying and amplifying the resultant signal, and using this 
signal to allow or block a fixed current from reaching an elec- 
tromagnet (Figure SI). With this sort of controlled magnetic 
field, gyrators become three terminal devices that act like gen- 
eralized transistors. Signals propagate from one side of the 
gyrator to the other (source to drain), but the gating current on 
the magnet controls the outgoing polarization. Importantly, 
treating the gate as the input and the drain as the output, the 
signal is amplified (as in electronic transistors). This is an 
advantage over isolator-based designs (e.g. some forms of op- 
tical computing [22]), as losses accumulate and hinder detec- 
tion. 

With these transistors, we can construct basic logic gates. 
First, consider a NOT gate. Given our basis, the easiest im- 
plementation is a 7r/2 gyrator (0 in gives 1 out, 1 in gives —0 
out. To remove the negative sign, which is an overall phase, 
the output can be fed into the source and gate of a tt gyrator, 
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FIG. 2: Implementing transistor logic with phonon gyrators. Ele- 
ments are denned in sidebar. X and Y are inputs, Z is output, and 
is an input fixed to logical 0. (a) NOT gate (b) NAND gate (c) OR 
gate (d) coherent XOR gate (using OR). 



the gyrator activated by 0). This design is good for chang- 
ing OR/NOR and AND/NAND, but an amplifying design is 
straightforward (Figure 2(a)). Logical is fed into the source 
of a gyrator and the signal to be inverted is applied at the gate. 
Logical at the gate activates the gyrator, rotating the source 
signal by ir/2 , yielding logical 1 at the drain. However, logi- 
cal 1 at the gate does not trigger the control, so the gyrator is 
not activated and the drain is at logical 0. Now consider two 
signals at the gate (Figure 2(b)). When the phonons are in a 
definite polarization, two signals will superimpose (no polar- 
ization change for identical inputs and ±7r/4 polarization for 
orthogonal signals). Control activates when this total signal 
has a horizontally polarized component greater than a critical 
value. Thus, if either input is logical 0, then control will acti- 
vate and the output will be logical 1. Only when both inputs 
are vertical will the output be logical 0. Therefore, this circuit 
is NAND. Control could alternatively select logical 1 (Figure 
2(c)). Now, the output will only be logical when both inputs 
are as well, ergo this circuit is OR. With either NAND/OR and 
NOT, we can use De Morgan's laws to create any logic gate. 
But for coherent phonons, XOR is simple (Figure 2(d)). Es- 
sentially, there is now destructive interference in the OR gate 
for signals of identical polarization, so the current can only 
exceed the critical amplitude when the signals are orthogonal. 
And since a ±7r/4 polarized signal activates control, this OR 
gate is now an XOR. 

While these designs are quite simple, there are other real- 
izations of logic gates. Exploiting De Morgan's laws, we can 
construct any operation using only 7r/2 gyrations and vertical 
control elements (Figure S2). Using gates in fixed arrange- 
ment and selectively (de)activating them, we can perform ar- 
bitrary operations (as in an FPGA). From this, single circuit 
parallel computing (RP) follows. Consider multiple signals of 
different frequencies all entering the same generic circuit. In- 
dependently controlling which gyrators and control elements 
work at which frequency, we can independently select the log- 
ical operation for each frequency. Such control is feasible 
with multistage gyrators (Figure 3) - several gyrators in se- 
ries where the sum of the gyrations at a given frequency is 
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FIG. 3: (COLOR ONLINE) Multistage gyrator: a concatenation of 
several gyrators. The lengths and magnetic field strengths at each 
stage are tuned to control the gyrations at multiple frequencies. 
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-B(j) is the magnetic field strength of the j stage, L(j\ its 
length, the frequency of the z t/l input, ^(oju\) it gyration, 
and <j> the gyration per unit length. The number of stages re- 
quired for a multistage gyrator scales as 0(N) for ungated 
elements and 0(N a ) (a G [1, 2] and depends upon the design 
used) for gated elements. These scalings can be improved by 
combining stages (allowing fringing fields to come to bear, 
slightly complicating the design and discussed in the supple- 
mentary materials). Control of different frequencies can be 
accomplished by filtering signals before amplification. 

Assuming uniform magnetic field, the polarization shifts 
obey Eq.l, one can easily find L^({B}, {<&}, {w}) by in- 
verting (f>(uju\, B(j\)({} denote sets). The exception to this is 
when the (f>(w(i\, B^) is separable ( <j> = f(oj)g(B)). While 
this is not generically the case with phonons, it is for optics 
(where <f> — V(ui)B) [23], the exceptions principally due to 
absorption (often in cold gases [24]). RP is (to our knowledge) 
presently unique to phonons (voltage is scalar, so manipulat- 
ing polarization is impossible in classical digital circuits), at 
least at room temperature. To show the feasibility of this type 
of design for phonons, we calculate Ln\{{B}, {$}, {w}) for 
a 2-frequency, 2-stage gyrator with no gyration induced at the 
lower frequency and a 7i72 gyration induced at the higher fre- 
quency. The gyration $(w; {B}, {L}) for the multistage gy- 
rator is plotted in Figure 4. Im[<i>]<0 denotes absorbance. 
Waveforms at the two frequencies are found analytically in 



real space (Figure 4 insets). The higher frequency mode has 
almost no loss (< 1CP 6 %), while the lower has a loss of 
19.8%. Finally, consider the number of stages possible for 
a given frequency and magnetic field strength range. We cal- 
culate (see supplementary materials) the maximum number 
of channels currently accessible is about 2700. Such a large 
number of channels is principally due to the large factor of 



Gyration vs Frequency for T uned Two-Stage Gyrator 





«> In] 



FIG. 4: (COLOR ONLINE) Rainbow Parallelization: designing a 2- 
stage gyrator to work at two frequencies. Total gyration as a function 
of frequency plotted in blue. Absorbtion (Im[$]) in dashed green. 
Circled frequencies have $ = and $ = nil. uj is rescaled by 
■q (the MA coupling rate), (inset) Real space representation of the 
waveforms for the two selected frequencies selected. L is rescaled by 
Cj P '/»7 (Cj P ' is the transverse phonon speed). Bottom/top figure is for 
lower/higher frequency mode with blue (orange) denoting horizontal 
(vertical) polarization. Vertical lines denote boundaries (2nd and 4th 
intervals are gyrators). 



YIG, which makes it possible to distinguish very small fre- 
quency separations. 

We have demonstrated a new mode of operation of an 
acoustic gyrator and used it to create basic logic gates. By ex- 
ploiting the wave nature of phonons and the tunable nonlinear- 
ity of MA, we demonstrate that MA materials can perform RP, 
an effect that is presently unique to this computing architec- 
ture. This suggests that there might be circumstances where a 
phonon computer could outperform an electronic computer. 

This material is based upon work supported by the Na- 
tional Science Foundation Graduate Research Fellowship un- 
der Grant No. 1 122374. Supporting data and figures relevant 
to this work are presented online under supplementary mate- 
rial. 
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I. CLASSES OF MAGNETO-ACOUSTIC MATERIALS 

There are several criteria for a material to be able to sup- 
port the magnetoacoustic interactions necessary for the acous- 
tic Faraday effect [13]. First, it must be a crystalline material 
with at least C3 (3-fold rotation) symmetry. Second, there 
must be some intermediate field to couple the phonon modes 
to the magnetic field. This is necessary because the Lorentz 
force acting on a phonon is generally negligible. The interme- 
diate field can be photons, electronic charge, or magnons. For 
photon-phonon coupling, the relevant interactions arise in low 
temperature metals, where a magnetic field allows the propa- 
gation of helicon (for uncompensated metals) and doppleron 
or Alfven wave (for both compensated and uncompensated 
metals) modes. These modes have been studied in pure metals 
(copper, nickel, etc), where it is found that dopplerons (which 
have phase and group velocities antiparallel to each other) are 
only supported for very small sets of wave vectors. 

In electron-phonon coupling, the relevant interaction is 
a coupling between the octupole moment of the 3d (as in 
CeA12), 4f (as in TGG), or higher order charge distributions 
and the crystalline electric field (CEF). These interactions 
arise because the nuclei in the crystal are charged, and so 
distortions of the lattice will induce distortions of the elec- 
tronic structure and vice versa. For materials of high rota- 
tional symmetry, the lower order modes of the charge distribu- 
tion (monopole, dipole, quadrupole) are less significant. Be- 
cause these materials have multiple phonon bands (i.e. acous- 
tic and optical branches), [15, 25] found that they have two 
types of magnetoacoustic interactions. In addition to the typ- 
ical coupling scheme (magnetic field to electronic charge to 
crystalline field to phonons) there is also an indirect coupling, 
where magnetoelastic effects in one branch can induce them 
in other branches due to phonon-phonon couplings. Addition- 
ally, there is a second class of electron-phonon interactions ca- 
pable of supporting magnetoacoustic effects. This is doppler- 
shifted cyclotron resonance (DSCR), which occurs under sim- 
ilar circumstances to helicon-phonon and doppleron-phonon 
couplings (i.e. low temperature pure metals). 

Lastly, in magnon-phonon coupling, the dependence of the 
exchange interactions in the spin Hamiltonian upon the lo- 
cation of the nuclei implies a dependence of the magnon 
bands upon the vibrational distortion of the lattice. This is 
essentially an ionic magnetostriction effect. It is found in 
both ferromagnets (for example YIG [12]), antiferromagnets 
or flopped antiferromagnets (as in Cr203 [14]), or paramag- 
nets doped with magnetic impurities (single molecule mag- 
nets, as in KMgF3+Ni [26]). It is an open question if phonon- 
phonon effects can induce a second type of magnetoacousic 



interaction for magnon-phonon coupling schemes (the con- 
clusions that we draw here still work qualitatively if these ef- 
fects are present, but the quantitative description would need 
to be modified). Since the acoustic Faraday effect is most 
extensively studied at room temperature in magnon-phonon 
coupling schemes (in particular in YIG), we shall YIG for our 
magnetoacoustic element. 

II. CONTROL MECHANISMS 

There are two basic classes of magnetic control. Control 
can either be direct or indirect. In direct control, the mag- 
netic field is modified by the phonon current itself, and in in- 
direct control the phonon current transduces some intermedi- 
ate signal which can control the magnetic field. Direct control 
can be performed by a ferromagnet (see Fig. Sl(a)). When 
a phonon current is incident upon the magnet, a portion can 
be absorbed or converted into localized modes that heat up 
the magnet. This increase in temperature will result in in- 
creased disorder of the magntization within the crystal, form- 
ing domains and weakening the total magnetization. Since the 
spontaneous magnetization is temperature dependent, the re- 
sultant magnetic field will be modified after a negligible delay 
(since changes in field propagate at the speed of light). If the 
temperature achieved is never so great as to demagnetize the 
entire magnet, there will be always be a net magnetization in 
a fixed direction, which means that upon cooling the down 
the direction of the total magnetization will remain fixed (thus 
keeping our gyrators acting like gyrators and not polarizers). 
While this approach is clearly the simplest to engineer, there 
are some limitations to direct control. First, if the thermal 
conductivity and specific heat of the magnet are not (approxi- 
mately) linear within the range of temperatures achieved, then 
it is exceedingly difficult to control the strength of the mag- 
netic field. What's more, the heating of the magnet depends 
upon the total thermal current, rather than any particular po- 
larization or frequency. While polarizers can eliminate trans- 
verse polarizations (and even frequeny intervals, when pairs of 
orthogonal polarizers are used), there is no easy way to elim- 
inate the longitudinal modes. While we can normally safely 
ignore these modes (as in every other element, they are com- 
pletely decoupled), here controlled operation would require 
us to either account for or eliminate any longitudinal effects. 
Lastly, because the phonons are being converted into magnetic 
disorder, they are being destroyed in this process and cannot 
be used in other circuit elements. 

Indirect control gets around many of these disadvantages. 
In the control scheme given here, there are three stages (Fig. 
S 1(b), stages shown in insets). First, the phonon current is sent 
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through a piezoelectric material (or some other linear, polar- 
ization sensitive transducer). The piezoelectric is cut such that 
the faces of the crystal correspond to the polarization vectors 
of the basis of transverse modes. The transverse modes will 
distort the crystal, inducing a fluctuating electrization which 
will produce an oscillating voltage gradient. This voltage gra- 
dient will be directed parallel to the polarization vector of the 
phonon current and have magnitude linearly dependent upon 
the phonon amplitude, that is d^V oc e iz (where e iz is the 
strain component in the ith direction of a mode propagating 
along z). Placing wires at the faces of piezoelectric will result 
in a very small AC voltage that depends solely upon a given 
phonon mode (i.e. AV oc fiie iz , where n is the normal vec- 
tor of the crystal surface). Because the energy extracted is so 
small, the phonon current leaves the piezoelectric essentially 
without loss. At the second stage, this voltage is converted 
into a finite, constant voltage (removing all non-transient time 
dependence is necessary to keep the magnetic field from fluc- 
tuating). This conversion is accomplished by using a concate- 
nation of voltage multiplying rectifiers. The number of recti- 
fiers used will increase the output voltage at this stage, so arbi- 
trary voltages can be achieved in principle. This voltage could 
be used to drive the current in an electromagnet (stage three), 
which would give a magnetic field that is approximately lin- 
ear in the amplitude of a polarized phonon current (linearity 
is maintained so long as the total magnetic field remains well 
below the maximum possible). Strict linearity, however, can 
be a drawback in building digital circuits. For example, in the 
construction of an OR gate it is preferable to have signals of 
strength e and 2e give identical fields (corresponding to inputs 
of X=l, Y=0 and X=l Y=l, respectively). To achive this sort 
of digital control, the voltage produced at the end of stage two 
is used to control the gate of an electronic transistor (see Fig. 
Sl(c)). Since the current produced is the one that flows from 
source to drain, the transistor output is (approximately) a dig- 
ital response to the gate voltage. The amplitude dependence 
of these three types of control regimes is summarized in Fig. 
Sl(d). Since large magnetic fields are typically required to 
achieve strong magnetoacoustic effects (typically in the range 
of .5 to 10 T [14, 21], sometimes fields of up to 20 T have been 
used [15]), the electromagnet will have to be quite powerful. 
Bitter electromagnets, for example, have been constructed to 
produce continuous fields of at least 35 T (superconducting 
magnets can exceed this, but not at room temperature) [27], 
far higher than any the fields required by our system. 
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FIG SI: Magnetic control designs, (a) Direct control by 
phonons. (b) Linear indirect control, (c) Digital indirect con- 
trol, (d) Magnetic field as a function of amplitude. The dotted 
line denotes direct control. The dashed line denotes linear in- 
direct control. The solid line denotes digital indirect control. 



III. SIMILARITY OF LOGIC GATES 



To determine the minimum number of gates required to 
construct all the two-terminal logic functions, we specialize 
to only using 7r/2 gyrations and veritcal control. By explicit 
construction (see Fig. S2), we see that we require three con- 
trol elements (one for each input, one for their superposition), 
and seven gyrators (three controlled, four fixed). Given that 
control elements are always paired to gyrators, this gives a to- 
tal of seven elements to perform all sixteen logic operations. 
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FIG S2: Alternative transistor logic using only vertical 
control and 90 degree gyrations. Implementations of all 16 
two-terminal logic functions. 



IV. INCORPORATION OF INHOMOGENEITY INTO 
MULTI-STAGE GYRATOR DESIGN 



The Faraday rotation induced by a uniform magnetic field is 
straightforward to calculate. When the magnetoacoustic inter- 
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action is due to magnon-phonon resonance, it takes the form 



/■ uj 2 

= L<j> 



n 2 



where $ is the total rotation, V is a constant (the analog 
of the Verdet constant), uj the phonon frequency, VL the ef- 
fective magnon frequency, z the runs over the thickness of 
the gyrator along the direction of propagation (L), and cf> is 
the gyration per unit length [13]. In a strong, uniform mag- 
netic field (B = Bz), Ml w -)\B\, where 7 is the gyro- 
magnetic ratio (g s /J,B for spins). When the magnetic field is 
not uniform (which can occur for multistage gyrators), calcu- 
lating the Faraday modes is more complicated. We assume 
a gradual inhomogeneity (that is, only including the effects 
of varying the magnon modes themselves, ignoring the de- 
viation of the coupled mode Hamiltonian from the Faraday 
regime). Consider the inhomogeneous magnetic field B = 
-\B'{z){x + y) + B{z)z, where B'{z) = d ^ L < B{z) 
and the term proportional to the derivative of the field is nec- 
essary to satisfy V ■ B = 0. This gives a vector potential 
A = —^B(z)(yx — xy), thereby maintaining the orbital Zee- 
man term proportional to B(z)L z . The spin Zeeman term, 
which is more important in this context, is proportional to 
B ■ S. The electronic degrees of freedom will arrange them- 
selves properly so that B(z)S z can be evaluated unambigu- 



ously at fixed z, i.e. that cj 
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is the magnon frequency for the nth branch/polarization and 
wavevector q, calculated to 0th order in transverse compo- 
nent) which can be integrated over in determining the Fara- 
day rotation. The terms perpendicular to z, on the other hand, 
will perturb the magnon Hamiltonian, shifting the frequen- 
cies. We treat these additional terms as a pertubation upon 

^nq (z), keeping terms less than 0{o?), as we did in our pre- 
vious derivation. The perturbation Hamiltonian, in Holstein- 
Primakoff coordinates is: 

U 



where 



X q B'(z)V2S{x-(-q)a(q) + x+(-q)a\q)} 



x ± (q) = -i[ 




F[rY^\f) 



{F is the Fourier transform operator, r is distance, and Y™ 
are the spherical harmonics), which (to lowest order) can be 
considered a constant so far as spin degrees of freedom are 
concerned. The first order correction vanishes trivially, and 
the second order is 

which is generically nonzero. 



V. DESIGN OF A MULTISTAGE GYRATOR 

While the use of inhomogeneous magnetic fields makes for 
a more elegant form of single circuit parallelization (in princi- 
ple, only a single magnetoacoustic material is required, with 
the magnetic field programmable by a set of control elements), 
it does make design more difficult. For uniform fields, the re- 
lation 

=E £cj)0( w (O' S Cj)) 
U) 

is clearly satisfied, meaning that we can invert the <f) matrix 
to determine the necessary length. To demonstrate how this 
design works, we consider the problem of designing a two 
frequency circuit which flips one mode and passes the other. 

fc 2 (l + 



To calculate the exact Faraday rotation we use k± 



uj ±L-7) ) wnere c| p 'fco = u> (cf is the transverse speed of 
sound), ui s the magnon frequency, and hr\ = 4|G '44 1 2 S 3 / 'M is 
magnetoelastic coupling energy (G44 is the coupling constant, 
S is spin, and M is the effective mass). The Faraday rotation 
equation given previously is based upon the assumption that 
fc+ - fc_ = {k\ - k 2 _)/2k (i.e. k+ + k_ w 2fc ). While this 
approximation is generally useful, it does not give an easy way 
to calculate the absorption of the incident wave. 

We work in units such that the equations simplify, i.e. 
il = lo s — i] w -fB — 77 + «r = B + iT, and rescale 
all the variables (w ->■ uj/rj , B -> B/r], T -> r/77, 



» 



k 2 



(p) 



giving the relation 



v± = u) (1 + B+ [ r±LU ) For clarity, we keep the magnon life- 
time (r s = 1/T) long, setting F =0.0023. We select Si=1.32 
and i?2=0.17, thus making it unlikely that our values happen 
to simplify the problem (we also want the magnetic field val- 
ues reasonably well separated, as significant losses appear in 
the range Re[cj] = [B, B + 1]). For $ 2 = 7r/2, it is best to 
select a frequency far from resonance, so we set ix>i=3.05. To 
make $1 = 0, it is necessary for it to be close to resonance 
with a magnon mode (this is implied by the shape of the Fara- 
day rotation curve). Thus, co 2 — 1-3 2 — T 2 . Solving for the 
lengths then gives L\ = 0.403 and L2=2.51. A plot of the 
total gyration as a function of frequency quickly confirms that 
this device acts as predicted. To translate this plot into the real 
space wavefunction, we assume impedance matching at the 
boundaries (i.e. ignoring losses due to reflection). This allows 
us to isolate the losses due to the finite lifetime, which we find 
tobe3.81x!0" 7 % atwi and 19.8% atw 2 - 



VI. THE MAXIMUM NUMBER OF CHANNELS 
CURRENTLY ACCESSIBLE 

The number of possible channels for a multi-stage gyrator 
is difficult to determine. In principle, an infinite number of 
stages could support an infinite number of channels, but since 
the complexity of the device increases as the number of stages 
increases. Moreover, there are limitations to the frequencies 
and field strengths that are feasible. As such, only a small sub- 
set of the channels physically accessible. Since the strength of 



4 



the magnetic field can be controlled experimentally so long as 
the electromagnets used are below their saturation magnetiza- 
tion, we are limited more by the resolution of phonon modes. 
We estimate the spectral resolution of the phonon modes us- 
ing the Q-factor. To make sure that the frequencies are dis- 
tinguisable, we assume ujn + = i^n+i — \^n+i 

where ujn is the Nth frequency used and Auj is the band- 
width. Essentially, this relation implies that frequencies are 
separated by their average bandwidth. Plugging in the defini- 
tion of Qjv = wat/Awn, we find that 

ui /uj = H N — - < ( ^ min iLl ^ 

mm x 

Therefore 

AT > ln( — )/ln(— -) « Q mi „ln — 

(the error in this approximation is <1% for Q>10). Setting 
cjjv = w max ,cJo = w min gives the number of distinguishable 



frequencies accessible. Experimentally, the majority of MA 
experiments concentrate in the range lj/2it e 10MHz to GHz 
[13]. For magnetic fields of <10T, however, the cyclotron fre- 
quency (i.e. the magnon frequency) is 152 MHz. We therefore 
lack the ability to easily tune gyrations above this frequency. 
Q for YIG is typically quite large in this frequency, values 
between 10 4 and 10 7 are regularly found for frequencies be- 
tween 10MHz and 1GHz (Q oc l/u) [28]. However, when the 
intensity of the phonon current is too great (lmW is the rule 
of thumb in the 500MHz-lGHz range), the magnon coupling 
is no longer linear and the Q factor degrades (the phonons lose 
energy to magnons, which quickly decay). This degradation 
is typically of order unity, but we shall take a conservative 
estimate of Q as 10 3 . This gives N =2700. This could be 
improved even further by countering the reduction of the Q 
factor, as well as extending the range of accessible frequen- 
cies. Conversely, given that this many channels would make 
fabrication dramatically more complicated, the a portion of 
the channels could be sacrificed (e.g. reducing the maximum 
strength magnetic field to a more readily achievable value). 
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